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Abstract: For high resolution imaging of a non-cooperative moving target, this 
paper proposes a sparse fusion imaging method. The imaging system contains two 
radar stations, which are separated by a certain bistatic angle and configured in a 
transmitter/receiver-receiver (T/R-R) manner. Consequently, two synthetic apertures are 
obtained at the same time from different aspect angles. By coherently fusing the echoes 
of the two radars, a virtual aperture spanned by these two sub-apertures can be constructed, 
which is larger than either of the sub-apertures; thus, the cross-range resolution of the image 
is enhanced. Moreover, the fusion of the echoes is realized by exploiting the sparse scattering 
property of the target. Then, based on the maximum a posteriori (MAP) criterion, the T/R-R 
fusion imaging problem is converted into a sparse signal recovery problem with unknown 
parameters. Finally, it is solved in an iterative manner, which contains two steps, i.e., sparse 
imaging and parameter estimation. Simulation results show that the proposed sparse fusion 
imaging method can improve the cross-range resolution significantly compared to inverse 
synthetic aperture radar (ISAR) within the same coherent processing interval (CPI). 

Keywords: sparse imaging; fusion imaging; T/R-R configuration; ISAR; bistatic ISAR; 
resolution improvement; moving target imaging 
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1. Introduction 

Microwave radar imaging is one of the major techniques for non-cooperative moving target 
recognition (NCTR). Compared to optical or infrared imaging, it has the ability to work under all 
weather, all time and long-range conditions. Therefore, it plays an important role in both defense and 
civilian applications [1-4]. 

Inverse synthetic aperture radar (ISAR) is the traditional microwave imaging method for a moving 
target [5-8]. In ISAR imaging, the target is mapped onto a slant range and cross-range plane, where 
the slant range resolution is determined by the bandwidth of the transmitted signal and the cross-range 
resolution is obtained by exploiting the motion of the target. Usually, the motion is decomposed into a 
translational and a rotational component [5]. After the translational motion is compensated, the rotational 
motion forms a synthetic aperture with size being A6 = tuT, where to is the rotational angular speed and 
T is the coherent processing interval (CPI). Since the cross-range resolution is inversely proportional 
to A9 and to is determined by the non-cooperative motion of the target, a long CPI time is required in 
order to obtain high-cross resolution. However, a long CPI time will make the motion compensation 
difficult [5-8]. As a variant of monostatic ISAR, bistatic ISAR (B-ISAR) has also been studied 
recently [9-12], where the transmitter and the receiver are separated by a bistatic angle. Although it 
can overcome some geometry limitations existing in monostatic ISAR, high cross-range resolution also 
cannot be guaranteed, due to the same reasons in ISAR. 



Figure 1. Imaging configuration. 




For the purpose of high resolution imaging with short CPI time, this paper proposes a sparse fusion 
imaging method by combining the monostatic ISAR and B-ISAR together. As shown in Figure 1, the 
imaging system consists of only two radar stations, which are separated by a certain bistatic angle 
and configured in a transmitter/receiver-receiver (T/R-R) manner. Consequently, two sets of echoes 
are received by Radar 1 and Radar 2, each of which corresponds to a sub-aperture located at different 
aspect angles. The interval between the two sub-apertures is equal to half the bistatic angle [9,10]. Under 
the condition of short CPI, usually, the half bistatic angle is larger than the sizes of the sub-apertures; in 
other words, there is an aperture gap between them. Therefore, by coherently fusing the received data, a 
virtual aperture spanned by these two sub-apertures can be constructed, which is much larger than either 
of the sub-apertures. As a result, the cross-range resolution will be improved significantly. 

In the echoes fusion procedure, due to relatively large aperture gap and short CPI, traditional 
Fourier-based sub-aperture fusion methods will suffer from high side lobes, and the spectrum 
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estimation-based gaped-data fusion methods [13,14] will also get performance degradation. In addition, 
three parameters (the rotational speeds of target relative to Radar 1 and Radar 2 and the bistatic angle) 
are required in order to fuse the echoes coherently. However, it is difficult to obtain accurate estimations 
directly from the return signals, due to short CPI time [15-17]. Here, we propose a new fusion imaging 
method by exploiting the sparse scattering property of the target. Based on the maximum a posteriori 
(MAP) criterion, the data fusion and the parameter estimation are combined together and converted into 
a joint optimization problem. Then, it is solved in an iterative manner by alternating two steps: sparse 
imaging and parameter estimation, where the first step is a sparse signal recovery (SSR) problem [18-22] 
and the second step can be solved by a linear search method. 

The remainder of the paper is organized as follows: Section 2 gives the imaging geometry and the 
return signal models. The basic idea of the imaging system and the sparse fusion imaging method are 
presented in Section 3. Then, the imaging algorithm is described in Section 4, and the numerical results 
are shown in Section 5 to validate the proposed method. Finally, Section 6 is the conclusion. 

2. Return Signal Modeling 

The imaging geometry is shown in Figure 2. Radar 1 transmits the signal and then receives the 
reflected signal, while Radar 2 works only as a receiver. For simplicity, we assume that the trajectory 
of the target is on the plane spanned by the line of sights (LOSs) of Radar 1 and Radar 2 during the 
short CPI time. It is worth noticing that the return signal model and imaging method introduced in 
this paper can be easily generalized to the case where the motion of the target is not parallel to that 
plane [23,24]. Let t m e [—T m /2, T m /2] be the slow time (snapshot time) and T m be the CPI. Figure 2 
shows two instants of the target at t m = 0 and t m > 0, where the target is represented with a dashed line 
and a solid line, respectively, and Ri(t m ) and R 2 (t m ) are the LOSs of Radar 1 and Radar 2. 

Figure 2. Imaging geometry. 

■ — Transmitted Signal 

y vA 
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radar 1 radar 2 

Tx-Rx Rx 



The space coordinate xOy is established with respect to the target by fixing the coordinate origin on 
the center of the target and setting Ri(0) (the LOS of Radar 1 at the snapshot time t m = 0) as the y-axis. 
According to the manner of ISAR, we define the y-axis as the slant range direction and the x-axis as the 
cross-range direction. 
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Let 0i(t m ) and f3 2 (t m ) be the angles of the two radars' LOSs with respect to the y-axis. Referring to 
our previous study [9], we treat the motions of the target relative to Radar 1 and Radar 2 separately and 
decompose each of the relative motions into a translational and a rotational component. During the short 
CPI time, the rotational motions of the target relative to Radar 1 and Radar 2 are approximately uniform 
rotations. Suppose the two rotational speeds are u>i and uj 2 , respectively; then (3\{t m ) and (3 2 (t m ) can be 
expressed as follows: 

Pl(t m ) = Ultra, foitm) = A)2 ~ U 2 t m (1) 

where /3 02 is the abbreviation of /3 2 (0), which is the bi static angle of Radar 2 at t m = 0. 

First we give the return signal model of Radar 2, because that of Radar 1 is the special case of Radar 
2. According to the high-frequency scattering mechanism, the target can be modeled as a set of dominant 
point scatterers [25-30]. Therefore, for the i th point scatterer located at = (x^, yi), under the condition 
of far field approximation, we can obtain the delay time at snapshot time t m as Equation (2): 

T2{t m ;xi,yi) « -[|Ri(t m )| + |R 2 (t m )| +yi-x i wa.p m , 
c 

+yiCOSf3 02 + XiUit m + XiU 2 tmCOSf3o2 + 2/iW2t m sin/3 0 2] (2) 

According to linear scattering theory, the return signal of Radar 2 is the superposition of the signals 
reflected from all of the scatterers. After matched filtering and translational motion compensation, the 
return signal of Radar 2 can be written as: 

sz{i, t m ) = 22 a i( x i, Vh i, t m )(Ji + e 2 (i, t m ) (3) 

i 

where t E [0, T p ] is the fast time, T p is the pulse repetition period, <7j is the scattering coefficient of the 
i th point scatterer, e 2 is the additive white Gaussian noise (AWGN) and a 2 is: 

\ 1 

t (Vi-Xi sin (3 02 + yi cos (3 02 ] 

c 



a 2 (xi, yu i, t m ) = sine InB 

x e _ j 27r c (Vi~ x i sm P02+Vi COS Po2+XiUJlt m +XiUJ2tm cos /3o2 +Vi^2t m sin ^02) 



where f 0 is the carrier frequency, B is the bandwidth and the migration through the range resolution cell 
is not considered. 

Then, we translate the return signal into its discrete form. Uniformly dividing the imaging scene 
into K x L spatial positions along x and y directions with grid sizes being dx and dy, respectively. As 
a result, the scattering coefficients of the imaging area can be expressed as a matrix 3 = \(Tk,i\KxL- 
Suppose that there are N samples in fast time and M samples in slow time; we use n to express the 
discretized fast time and m to represent the discretized snapshot time. Then, the return signal of Radar 2 
can be expressed as an N x M matrix: 



K L 

S 2 = ^^A^ct m +E 2 k = l,...,K, l = l,...,L 

k=l 1=1 



(5) 
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where A 2 ' = [a 2 ' n m ] jvxm is the observing matrix with respect to the point scatterer located at (k, I) and 
E 2 is the N x M noise matrix. 

Next, we express Equation (5) in its vector form: First let s 2 = vec(S2), where vec(-) refers to the 
vectorization operation (i.e., stacking the columns of a matrix on top of each other) [22]. Then, define 
the observing matrix of Radar 2 as: 



vec (A2' 1 ) , vec (A 2 ' 1 ) , vec (Af' 1 ) , vec (A2' 2 ) , vec (A2' 2 ) , vec ^Af' 2 J 



...,vec ( Ai' L ) ,vec (A 2 2 ' L ) , vec (A 



K,L 
L 2 



(6) 



Furthermore, let cr = vec (S) and e 2 = vec(E 2 ). In this way, the return signal of Radar 2 can be 
expressed in the following compact form: 



s 2 = A 2 er + e 2 



(7) 



The return signal of Radar 1 is the special case of that of Radar 2. Let (3 02 = 0 and replace cu 2 with 
Ui, we can obtain: 



Si = Aicr + ei 

where Ai is the observing matrix of Radar 1 and ei is an AWGN vector. 



(8) 



3. Sparse Fusion Imaging in T/R-R Configuration 

Figure 3 shows the basic idea of the proposed imaging method from the perspective of synthetic 
aperture, where A61 and A9 2 are the apertures of Radar 1 and Radar 2, respectively. In a T/R-R imaging 
system, the two radars receive return signals from different viewing directions, so A9i and A9 2 are 
obtained at the same time and located at different aspect angles. Suppose the interval between them is 
do; then, according to the theory of ISAR and B-ISAR [5,9-1 1], we can obtain that: 

A0i = T m ui 

AO 2 = T m (ui + u 2 cos ^02) /2 

d e = (3 02 /2 (9) 

Figure 3. Schematic diagram of the synthetic apertures of the transmitter/receiver-receiver 
(T/R-R) imaging system. 

aperture of radar 1 aperture gap aperture of radar 2 




virtual aperture after fusion 



As is well known, large synthetic aperture is required in order to improve the cross-range resolution 
of the image. In this paper, we will try to construct the large virtual aperture Ad by coherently fusing the 
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echoes Si and s 2 . As shown in Figure 3, this virtual aperture spanned by the sub-apertures A6i and A6> 2 
has a size of: 

A6 = (A6i + A9 2 ) /2 + d e (10) 

Because dg = (3 02 /2 > 0, AO is larger than either of A9\ and A6 2 . Specifically, we consider two cases: 

(1) when 0 < dg < (A6 1 ! + A8 2 ) /2, A6i and A9 2 will overlap with each other partially, and we 
can obtain that (A6>i + A9 2 ) /2 < AO < (A8i + A6 2 ). In this case, the cross-range resolution can be 
improved by fusing s 1 and s 2 , but the improvement factor will be less than two. 

(2) when d e > (A6>x + A9 2 ) /2, the two sub-apertures will be separated without overlapping; in 
other words, there will be an aperture gap between them, as shown in Figure 3. Especially under the 
condition of short CPI and a relatively large bistatic angle, the aperture gap will be much larger than the 
sub-apertures themselves. As a result, AO will become much larger than (A6*i + A9 2 ). If this virtual large 
aperture can be obtained, then the cross-range resolution will be improved significantly. Furthermore, 
note that in order to obtain AO by coherently fusing s 1 and s 2 , the scattering properties of the target with 
respect to the LOSs of Radar 1 and Radar 2 must be approximately the same. Therefore, generally, the 
bistatic angle (3 02 is constrained to be less than 10° [31]. 

For the first case, the spanned aperture can be obtained by traditional sub-aperture fusion 
methods [23,24]. This paper will focus on the second case, which is more likely to occur in T/R-R 
fusion imaging for non-cooperative moving target with short CPI. Above all, based on the return signal 
models (7) and (8), we combine them together and construct a fusion model as follows: 

s = Acr + e (11) 

where: 



A= ' , e= 1 (12) 






s is the combined return signal vector; A is the corresponding observing matrix, and e is the noise vector. 
Suppose ei and e 2 are independent and identically distributed AWGN vectors; then, e is also an AWGN 
vector with the mean being zero and the covariance matrix being rjl, where 77 is the power of the noise, 
which is an unknown parameter to be estimated. 

According to Equations (4) and (1 1), A is determined by three unknown parameters oji, u 2 and (3 02 . 
To fuse Si and s 2 coherently, these three parameters must be estimated. However, due to the short CPI 
time, it is difficult to obtain accurate estimations of the parameters directly from the return signals in 
advance [15-17]. In this paper, we combine the fusion imaging procedure, as well as parameter 
estimation together and treat the scattering coefficients vector cr also as an unknown parameter. Then, 
there are a total of five unknown parameters in Equation (11), i.e., cr, 77, loi, u 2 and (3q 2 . Therefore, 
the fusion imaging problem (11) can be seen as a multi-parameters estimation problem. Assuming that 
the parameters cr, 77, U\, co 2 and (3 02 are uncorrected, then according to the MAP criterion [32], we can 
obtain these parameters by the following joint optimization problem: 



(c>, 77, 0)1,0)2, ^02 J = arg max / (cr, 77, u lt u 2 , /3 0 2 |s 



(T,r),uii,u>2,Po2 (13) 

ocarg max / (s |cr, 77, u u u 2 , P02) f (cr) / (77) / (wi) / (w 2 ) / (Ah) 
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where cr, r), &\, Cj 2 , P02 are the estimations of cr, 77, w ls w 2 and /3 0 2, respectively. 

Because e is an AWGN vector, the conditional probability density function of s is: 

/(s|rr,77,Wi,a;2,A)2) = CAf (Acr , r]l) (14) 

where CM represents the complex Gaussian distribution. 

As there is no a priori knowledge about 77, cu 1 , cu 2 and (3 02 , thus we use the non-informative prior [33] 
as their distributions according to Bayesian statistics: 

/ iv) f(ui)f(U2)f(P02) OC 1, 77 E [0, f] max ],UJi E [0, Wi )max ], W 2 e [0, W 2 , max ], A)2 e [0, A)2,max] (15) 

where ?7 max , wi, max , W2,max and /3 0 2,max are the maximums of 77, wi, u; 2 and /3 02 , respectively. 

If we suppose cr also as a non-informative variable, Equation (13) will become a traditional minimum 
mean square error (MMSE) estimation problem [32]. In this case, the solution of cr will suffer from high 
side lobes and low resolution due to the relatively large aperture gap and short CPI [34], even if 77, u>i, 
co 2 and (3 0 2 are accurately known, not to mention that they are actually unknown parameters. In order 
to solve Equation (13) efficiently, here, we exploit the sparse scattering property of the target under the 
condition of high-frequency scattering, and then, assume the sparse a priori distribution of cr as [22]: 

K L 

/(<r)«nn e " f(|Bwl '~ i) (i6) 

fc=i ;=i 

where q E (0, 1] is a user parameter. When q — y 0, the peak of f(cr) at cr = 0 becomes infinite. This 
means that cr will be equal to zero with a high probability, namely cr will have a sparse solution in 
Bayesian inference. 

Substituting Equations (14)— (16) into Equation (13) and taking the negative logarithm form of 
Equation (13), we can obtain: 

(<r, 77,0)1,0)2,302 J = arg min C (cr, 77, u x , u 2 , Pw) 

V / CT,77,£ji,aJ2,/3o2 

= arg min { MN\og V + ± ||s - A (u h w 2 , P02) *\\ 2 2 + £ E f (M ~ 1) 

where C represents the cost function and ||-|| 2 denotes the Euclidean norm of a vector. 

Therefore, based on the MAP criterion and the sparse scattering property of the target, we have 
converted the T/R-R fusion imaging problem (11) into a joint optimization problem (17). The algorithm 
of Equation (17) will be discussed in the following section. 



4. Imaging Algorithm 

The purpose of Equation (17) is to find the sparse solution of cr; hence, it is an SSR problem about cr. 
In the traditional SSR problem, the observing matrix A usually is known in advance, and all of the efforts 
are focused on the recovery of the sparse solution [22]. However, in Equation (17), the observing matrix 
A contains three unknown parameters. According to [25,35], here, we propose a cyclic optimization 
method to solve Equation (17), which recovers the sparse solution, as well as estimates the parameters in 
an iterative manner. Firstly, we fix the three parameters, and then, Equation (17) becomes the usual SSR 
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problem, which can be solved by sparse learning via the iterative minimization (SLIM) method proposed 
in [22]. Secondly, we fix <x and 77, and estimate ui\, uj 2 and (3 02 by minimizing the cost function C. Then, 
repeat the two steps until the solution converges. The details are described as follows. 



4.1. Algorithm Description 

The initial values are set as follows: <x° is set as traditional range-Doppler (RD) imaging result [10]; 
uji, Cj 2 and (3q 2 are estimated roughly by imaging geometry and the motion of the target; and 
r)° = jjjj ||s — A 0 *?- 0 ^. We assume that fr\ fj l , Cj\, Cj % 2 and (3 l Q2 are the estimated values after the i th 
iteration; then, the (i + \) th iteration is shown as follows: 



Step 1 : Sparse Imaging [22] : 

The purpose of this step is to estimate <x l+1 and fj l+1 while fixing Cj\, Cj\ and /?q 2 . By using these 
three fixed parameters, the observing matrix A 1 (uj\, co 2 , (3 l 02 ) can be determined. Then, the optimization 
problem (17) becomes: 

(& i+1 ,f> i+1 ) = argminCfo-,//) (18) 

which has the same form as the usual SSR problem. It can be solved by the SLIM algorithm [22], which 
is also a cyclic optimization procedure. When the SLIM algorithm converges, we obtain <x i+1 and fj l+1 . 



Step 2: Parameter Estimation: 



In this step, we try to estimate u 2 +1 and /^J 1 by fixing <x i+1 and fj l+1 . Removing the terms 
which have nothing to do with these three parameters, then the optimization problem (17) becomes: 

(ui +1 ,u t 2 +1 J t 0 t 1 )=aTg min C p (u u u 2 , Am) 

= arg mm s - A (u) 1 , u 2 , p 02 ) a 

where A (u>i, ui 2 , /3q 2 ) cr' L+1 is the estimation of the return signal based on <r 4+1 . The cost function C p is 
the error between A (ui, u 2 , /3 02 ) and s. According to our study, C p u 2 , /3 02 ) is a non-convex 
function, so we propose a three-dimensional linear search method to solve Equation (19). 
First, we define the discrete search space as A: 

A = {(wi, cu 2 , (3 02 ) \ui e lin Ai) , u 2 G lin (fi 2 , A 2 ) , A)2 e Hn (fi 02 , A 02 )} (20) 

where fi 2 an d £^02 represent the search range of U\, co 2 and (3 02 , A 1; A 2 and A 02 are their 
corresponding step sizes, function lin (£2, A) generates a linear discrete space, whose range is £2 and 
step size is A. Then, by searching the minimum of C v (ui\, co 2 , (3 02 ) on A, we will obtain the optimal 
estimations of these three parameters. 
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Step 3: Convergence Judgment: 

After Step 1 and Step 2, we can obtain cr %+1 , r) l+1 , u\ +l , Cj 1 ^ 1 and /3q2 1 - Then, iteratively operate the 
two steps until reaching the stop criterion: 

||^ +1 -^|| 2 /||^ +1 || 2 <7 CT (2D 
where 7 CT is a small positive number. 

4.2. Some Discussions about the Algorithm 

In this subsection, we will give some discussions in terms of fast implementation, computational 
complexity and convergence. 

(1) Methods to decrease the computational burden: Because <x is a sparse solution, this means that 
most elements of er are zero or close to zero. Therefore, in Step 2, we only select the large elements of a 
to calculate the cost function C. In this way, the computational burden of Equation (19) will be reduced 
significantly. 

(2) Computational complexity: According to [22] the major computation of Step 1 comes from the 
matrix vector product (MVP). Obviously, in Step 2, the major computation of Equation (19) is also 
contributed by the MVP. Therefore, MVP is the main computation of the whole algorithm, and the 
computational complexity is O(MNKL). 

(3) Convergence of the algorithm: In order to guarantee the convergence of the proposed algorithm, in 
this paper, we adopt two measures. First, the SLIM algorithm [22] is exploited in Step 1. SLIM belongs 
to the Bayesian compressive sensing (CS) method, which is relatively robust against noise. Second, in 
Step 2, only the strong scatterers are selected to calculate the cost function C, which can reduce the 
impact of noise on the convergence of the algorithm to a certain degree. Although, it is difficult to prove 
the convergence of the algorithm directly. However, by the above two measures, we can ensure that the 
cost function will decrease after every iteration. Moreover, according to the simulation results, usually 
the cost function will converge after only a few iterations. 

5. Simulation Results 

In this section, we demonstrate some simulation results to verify the effectiveness of the proposed 
T/R-R fusion imaging method. The simulation conditions are as follows: the target shown in Figure 4 
is modeled as a set of dominant point scatterers with different scattering coefficients. During the short 
CPI, it is moving along the x-axis at the speed of v = 300 m/s. Assume that the bistatic angle is 
[3o2 = 4.3° at t m = 0 and the distances between the target and Radar 1 and Radar 2 are 
Ri (0) = 10 km and R 2 (0) = 12 km, respectively. Then, according to the theory of ISAR and 
B-ISAR, we can calculate that the truth-values of the rotational speeds are lui = 0.030 rad/s and 
co>2 = 0.025 rad/s. The signal-to-noise ratio (SNR) of the return signal is 15 dB. More details about 
the simulation conditions are given in Table 1, and the user parameter q in Equation (16) is chosen as 0.6 
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(in fact, we found not much of a difference of the various values of q < 1). In the simulation, the initial 
values of ui, u 2 and /3 02 are set to 0.04 rad/s, 0.04 rad/s and 6°, respectively. 

Before imaging simulation, we calculate the theoretical resolutions of Radar 1, Radar 2 and the 
T/R-R fusion system according to the theory of ISAR and B-ISAR. In the slant range direction, the 
resolutions are determined by the bandwidth of the transmitted signal and the bistatic angle, as shown in 
Table 2. Because (3 0 2 = 4.3° is small, they have approximately the same slant range resolutions 
as 0.38 m. The cross-range resolutions are inversely proportional to the sizes of the apertures A6*i, 
A6*2 and A9. These apertures can be obtained by substituting the truth-values of cui, co 2 and [3o 2 into 
Equations (9) and (10), where T m = 0.25 s. Then, the cross-range resolutions are calculated and shown 
in Table 2, where A = 0.03 m is the wavelength corresponding to the carrier frequency. We can see that 
if Radar 1 and Radar 2 are considered separately, their cross-range resolutions are 2.00 m and 2.20 m, 
respectively, which are very poor compared to their slant range resolutions. The proposed T/R-R fusion 
imaging system constructed a large spanned aperture with the size being 2.56°, which is much larger 
than A#i and A9 2 . Thus, the cross-range resolution is improved to 0.34 m, which is more than five-times 
better than either Radar 1 or Radar 2. 

Figure 4. Target model. 




Cross Range (m) 



Table 1. Simulation Conditions. 



System Parameter 



Value 



Carrier frequency of the signal fo 
Bandwidth of the transmitted signal B 
Pulse repetition frequency PRF 
Numbers of samples of fast-time and slow-time 
Image size 

Separations between cross- and slant range bins 
Numbers of cross- and slant range bins 



10 GHz 
400 MHz 
80 Hz 
N = 30 M = 20 

Length = 12 m, Width = 6 m 
dx = 0.2 m, dy = 0.2 m 
K = 60, L = 30 
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Table 2. Aperture sizes and theoretical resolutions. 



Slant Range Resolution Aperture Size Cross-Range Resolution 



Radar 1 



T/R-R 



psr.l 



25 = 0.38 m 



Radar 2 p sr , 2 = 2Bco J {M2) ~ 0-38 m 



P. 



sr — 2gcos 2 (/3 02 /4) 



0.38 m 



A6»i = 0.43° 
A6 2 = 0.39° 
A6> = 2.56° 



Per, X 
Per, 2 
Per 



A 



2A0i 

A 
2A0 2 

A _ 
2A0 " 



= 2.00 m 
= 2.20 m 
0.34 m 



5.1. Imaging Simulation 

The proposed sparse fusion imaging method is utilized to reconstruct the image of the target. For 
the purpose of comparison, the images obtained by Radar 1 and Radar 2 separately are also given. 
Furthermore, all of these three kinds of images are reconstructed both by the traditional RD method and 
the SSR method. The imaging results are shown in Figure 5. 

Figure 5. Imaging results. In (a)-(e), the parameters uo\, tu 2 and (3 02 are set to be 
the truth-values in advance, because the range-Doppler (RD) method or the sparse signal 
recovery (SSR) method with the short coherent processing interval (CPI) does not have the 
ability to estimate the parameters accurately, while in (f), the parameters are estimated during 
the imaging iteration by our method. 
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(a)Image of Radar 1 by the RD method. (b)Image of Radar 2 by the RD method. (c)Image of T/R-R by the RD method. 
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(d)Image of Radar 1 by the SSR (e)Image of Radar 2 by the SSR (f)Image of the proposed method, 
method. method. 



Figure 5a shows the image of Radar 1 by the RD method. One can see that in slant range direction, 
it has the ability to distinguish the point scatterers, because of the sufficient signal bandwidth. However, 
in cross-range direction, its resolution is too poor to distinguish the point scatterers. Figure 5d is 
the SSR result corresponding to Figure 5a. In this figure, only a small part of the point scatterers 
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are recovered with accurate locations and amplitudes; many scatterers are recovered with errors in 
locations or amplitudes, and the others are missed. That is because the aperture of Radar 1 is 
so small that even the SSR method cannot recover the image correctly. The images of Radar 2 
obtained by the RD and SSR method are shown as Figure 5b,e, respectively. Similar to Radar 1, 
neither of the methods can reconstruct the image correctly, due to the small aperture. Figure 5c 
shows the image obtained by directly solving the fusion model (11) with the RD method, where 
the parameters are set to be their truth-values. It shows that although the cross-range resolution 
is enough to distinguish the point scatterers; the image suffers from very high sidelobes, due to 
the large aperture gap between the two sub-apertures. Finally, the image obtained by the proposed 
sparse fusion imaging method is shown in Figure 5f. Compared to the true target model shown in 
Figure 4, we can see that the image is well focused with high resolution and low side lobes; both the 
locations and amplitudes of all the scatterers are recovered correctly. Therefore, it proves that our method 
has the ability to reconstruct a high resolution image with a very short CPI time. 

5.2. Algorithm Performance Simulation 

The convergence, as well as the parameter estimation performance of the algorithm are assessed 
by 60 Monte-Carlo trials. Figure 6a shows the convergence performance of the cost function under 
different SNR conditions. We can see that when SNR = 0 dB, the cost function cannot converge, due 
to the strong noise, but when SNR > 5 dB, the cost function will converge after only a few iterations. 
Moreover, with the increase of the SNR, the convergence rate becomes faster. In order to represent the 
recovery performance of the proposed method under different SNR conditions, we define the normalized 
root mean square error (NRMSE) between the imaging result and the original target model as follows: 

NRMSE =||cr m — cr|| 2 /||cr|| 2 (22) 

where ||-|| 2 denotes the Euclidean norm, cr is the scattering coefficient vector of the target model and 
cr m is the average estimation of cr over 60 Monte-Carlo trials. NRMSE demonstrates the recovery 
performance of the imaging method; a smaller NRMSE means better recovery performance. Figure 6b 
shows the values of NRMSE versus SNR. It shows that with the increase of the SNR, the values of 
NRMSE decrease, i.e., the recovery performance becomes better. When SNR > 5 dB, the NRMSE is 
close to zero, which means that under this condition, the imaging result is substantially the same as the 
original target model. 

Figure 7 shows the parameter estimation accuracy, Figure 7a,b is the means of ui, u 2 and /3 02 versus 
SNR, and Figure 7c,d is their variances. It shows that when SNR > 5 dB, the average estimations of u>i, 
C02 and (3q2 are very close to their truth-values: 0.030 rad/s, 0.025 rad/s and 4.3°, respectively, and the 
variances of the estimations are fairly small. This indicates that the proposed method has high parameter 
estimation accuracy. 
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Figure 6. (a) Convergence performance of the algorithm under different SNR; 
(b) normalized root mean square error (NRMSE) of the images versus SNR. 
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6. Conclusions 

This paper has proposed a sparsity-based fusing imaging method for a moving target in T/R-R 
configuration. Based on the idea of coherently fusing two widely separated apertures into a large aperture 
and by exploiting the sparse scattering property of the target, the cross-range resolution of the image is 
improved significantly. The details of the method, as well as its corresponding algorithm are presented. 
We also give some comments on the fast implementation and the convergence of the algorithm. Finally, 
some simulation results are given to validate the proposed method. They show that the proposed sparse 
fusion imaging method has the ability to reconstruct a high resolution image with a very short CPI time, 
and the algorithm can converge within a few iterations. 
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